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The recently reported positron fraction up to ~ 350 GeV by AMS-02 seems to have tension with 
the total electron/positron spectra detected by Fermi and HESS, for either pulsar or dark matter 
annihilation/decay scenario as the primary positron sources. In this work we will show that the 
tension will be removed by an adjustment of the primary electron spectrum. If the primary electron 
spectrum becomes harder above ~ 50 GeV, similar as the cosmic ray nuclei spectrum, the AMS- 
02 positron fraction and Fermi/HESS data can be well fitted by both the pulsar and dark matter 
models. This result indicates that there should be a common origin of the cosmic ray nuclei and 
the primary electrons. Furthermore, this study also implies that the properties of the extra sources 
derived from the fitting to the AMS-02 data should depend on the form of background. 
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I. INTRODUCTION 

The AMS-02 collaboration reported the very precise 
measurement of the positron fraction e + /e with ener- 
gies up to 350 GeV recently [l|. The positron fraction 
shows a continous increase up to ~ 100 GeV, which is 
consistent with the previous PAMELA result [2, Q and is 
lower than that measured with Fermi-LAT |4| ■ A flatten- 
ing of the positron fraction above ~ 100 GeV is revealed 
by the AMS-02 data, for the first time. The AMS-02 data 
implies that there is excess of positrons above tens of 
GeV compared with the standard cosmic ray (CR) back- 
ground, and the amount of excess positrons should be 
less than previously estimated according to the PAMELA 
data. 

Several works appears to explain the AMS-02 data 
with pulsars Q or dark matter (DM) scenarios |g-|l|. A 
thorough study of the properties of the extra positron 
sources, including the astrophysical one like pulsars and 
the dark matter (DM) scenario, based on the AMS- 
02 data and the electron (or e ) spectra measured by 
PAMELA @, Fermi-LAT (HlHU and HESS dHH, was 
given shortly after the publication of the AMS-02 data 
[14J (Paper I). It was found that there was difficulty to 
fit the AMS-02 positron data and the Fermi-LAT/HESS 
total electron spectra simultaneously, either in the pul- 
sar scenario or in the DM scenario. The results seem to 
imply that there might be tension between the AMS-02 
data and the Fermi-LAT/HESS data, in the present the- 
oretical framework. This conclusion has been confirmed 
by other studies [15|, [lj|- One possible reason leading 
to the tension is the constraint on the electron injection 
parameters by the pure electron spectrum by PAMELA. 
If the PAMELA data are not included, the primary elec- 
tron spectrum has larger free space and the AMS-02 data 
and Fermi data can be easier to be fitted simultaneously, 
as shown in some recent works to explain the AMS-02 



result 

Several possibilities to reconcile these two data sets 
were discussed in Paper I, including multiple components 
of the extra sources and the existence of spectral harden- 
ing of the primary electron spectrum. The basic idea of 
spectrum hardening is like this. If the primary electron 
spectrum becomes hard at high energy we need a hard 
electron/positron spectrum from the extra source to give 
correct positron fraction. This also leads to higher e 
spectrum to be consistent with Fermi data. The spec- 
trum hardening was stimulated by the observed spectral 
hardening of the nucl ei sp ectra in recent years by sev- 
eral collaborations |17H19| . A unified spectral hardening 
at rigidity R ~ 200 GV (or E k ~ 200 GeWn)was mea- 
sured precisely by PAMELA or CREAM [l8|,[ig|. If there 
is a spectral hardening of the CR nuclei spectra, it is nat- 
ural to expect a similar hardening of the primary electron 
spectrum. 

Models to explain the spectral hardening include the 
multi-component sources [20h24| . non-linear acceleration 
of the particles [25| , or the propagation effect [26|, (27[ . In 
[28f | the authors pointed out that if there was a spectral 
hardening of the primary electron spectrum, there would 
be a less steep increase (or decrease) of the positron frac- 
tion above - 200 GeV. 

In the work we investigate in detail whether to involve 
such a spectral hardening of the primary electron spec- 
trum can help eliminate the tension between the AMS- 
02 data and the Fermi-LAT/HESS data. We employ the 
CosRayMC tool developed in [2!| to fit the observational 
data within the high dimensional parameter space. The 
GALPROP package 1 [3f| which has been embedded in 
CosRayMC is used to calculate the propagation of the 
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charged CRs. The diffusion-reaccclcration propagation 
frame is adopted, and the major propagation parame- 
ters arc A)|fl =4GV = 5.94 x 10 28 cm 2 s~\ S = 0.377, 
va = 36.4 km s _1 and Zh = 4.04 kpc [3l|. The goodness 
of fit, constraints and implication of the model parame- 
ters are discussed. 

In the next Section we simply describe the models to 
fit the data. The results are presented in Sec. III. In Sec. 
IV we give the discussions and conclusions. 



II. MODEL 

In this section we describe the major aspects of the the- 
oretical models to reproduce the electron/positron data 
briefly. The injection spectra of the primary protons 
(heavier nuclei are less important in this study) and elec- 
trons are both assumed to be broken power-law functions 
with respect to momentum p 



q(p) <x <&>/<!)-», 



P < Pbr,l 
p,e 



<i<p<< 2 (i) 



>/< 2 )- T2 « 2 /<i)- T1 , P>vZa 

where Pbr,i represents the low energy break, pbr,2 is the 
high energy break to be responsible for the spectral hard- 
ening, 70, 71 and 72 are the spectral indices in different 
momentum ranges. We also employ the log-parabolic 
function to describe the spectral hardening of the elec- 
trons, i.e., q( P ) OC (p/p^J-Ti+Talogfr/MeV) for p > p ^^ 

In this case p^ v 2 is not used. The absolute fluxes of pro- 
tons and electrons are determined through normalizing 
the propagated fluxes to normalization factors A p and 
A e . 

The background positrons are expected to be produced 
through the collision of CR nuclei with the interstellar 
medium (ISM) during t he p ropagation. The parameteri- 
zation of pp collision in |32j is employed to calculate the 
secondary production of positrons and electrons. Similar 
as done in Paper I, we further introduce a free factor c e + 
to adjust the absolute fluxes of the secondary positrons 
and electrons to fit the data. Such a factor may repre- 
sent the uncertainties of the hadronic interactions, prop- 
agation models, the ISM density distributions, and the 
nuclear enhancement factor from heavy elements. 

In the PAMELA era it was found that the background 
contribution are not enough to ex plai n the observed 
positron fraction and total e^ data [33|, [34[ . Therefore 
the extra sources of e ± beyond the traditional CR back- 
ground are introduced to explain the data. We will base 
on the same theoretical framework to fit the AMS-02 data 
in the work, assuming continuously distributed pulsars or 
the DM annihilation/decay to be the extra sources of e . 

The injection spectrum of e^ from pulsars is assumed 
to be power-law with an exponential cutoff 



where A psr is the normalization factor, a is the spectral 
index and p c is the cutoff momentum. The spectral index 
a is limited in the range 1.4 to 2.2 according to the 7-ray 
observations of pulsars [35| . The spatial distribution of 
pulsars is taken to be the cylindrically symmetric form 
given in [36| 
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(3) 
where R@ = 8.5 kpc is the distance of the solar location 
to the Galactic center, z s « 0.2 kpc is the characteristic 
height of the Galactic disk. 

As for DM scenario (taking annihilation as illus- 
tration), we focus on the leptonic two-body annihila- 
tion channels /i + /i~ and t + t~, as implied according 
to the PAMELA and Fermi-L AT/HESS data of the 
electrons/positrons and the antiprotons [37H39| . The 
positron/electron production function from DM annihi- 
lation is (assumed to be Majorana particles) 



q{r,p) 



(av) dN 
2m 2 dp 



x p 2 (r), 



(4) 



where m x is the mass of DM particle, (av) is the velocity 
weighted annihilation cross section, dN / dp is the yield 
spectrum for one annihilation of a pair of DM particles, 
and p(r) is the DM density profile. The spatial profile of 
DM energy density is taken to be Navarro-Frenk- White 
(NFW, J4fl|) distribution 



p{r) 



{r/r s )(l+r/r s y 



(5) 



with parameters r s = 20 kpc and p s = 0.26 GeV cm~ 3 . 

For low energy particles we further employ a simple 
force field approximation to take into account the so- 
lar modulation effect [4l|. Since the operation period of 
PAMELA and AMS-02 is close to the solar minimum, 
the modulation potential is required to be smaller than 1 
GV. Note, however, the low energy part of the positron 
fraction measured by PAMELA and AMS-02 might not 
be reproduced with such single solar modulation model, 
and more complicated charge-sign dependent modulation 
effect is necessary (43, |43| . 



III. RESULTS 



We first determine the parameters of the proton in 



q(p) = A psl p Q cxp(-p/p c ), 



jection spectrum through fitting to the PAMELA [18 1 
and CREAM [3 data. For CREAM data we include 
10% systematic uncertainties as discussed in [I3. The 
high energy break p^ T 2 is fixed to be 230 GeV as sug- 
gested by the PAMELA data. The best fitting parame- 
ters of the proton spectrum are: 70 = 1.80, 71 = 2.42, 
72 = 2.33, p£ r 1 = 12.3 GeV, and solar modulation po- 
tential </> = 495 MV. The normalization of the prop- 
(2) agated proton flux at 100 GeV is A p = 4.55 x 10" 9 
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Comparison of the best fitting 



spectrum of protons with the observational data is shown 
in Fig. [TJ We see very good agreement between the cal- 
culated spectrum and the data. The minimum x 2 value 
is about 24 for 72 degree of freedom (dof). 
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FIG. 1: Proton spectrum derived through fitting PAMELA 
and CREAM data. References of the proton data: AMS [HI, 
BESS HJ|, ATIC2 0, PAMELA QJ] and CREAM QJ]. 



Since the observational period of protons by PAMELA 
is almost the same with that of electrons by PAMELA 
and positrons by AMS-02, we should expect a common 
modulation amplitude for these particles (besides the 
charge-sign dependent effect). Therefore we employ a 
prior on the modulation potential <j> = 500±53 MV comes 
from the fit of the proton data. 

We first fix the electron second break energy p§ r 2 at 
230 GeV, same as that of protons. The best fitting results 
of the positron fraction and electron spectrum are shown 
in Fig. [2] The fitting parameters and the x 2 value are 
presented in Table |U Compared with the case without 
spectral hardening of the primary electrons (Paper I), 
the fitting is indeed improved. The x 2 value decrease 
from ~ 280 to 246 with one more parameter. However, 
the overall fitting is still not satisfactory. We can see from 
Fig. [2] that when AMS-02 data are well reproduced, the 
model expectation is lower than the Fermi data, which is 
similar with the findings in Paper I. 

We then relax the break momentum of the electrons 
and redo the fit. In this case we find the improvement is 
significantly, as shown in Fig. [3l The parameters are also 
given in Table HI The minimum x 2 value over dof is about 
1.1, which implies a rather good fitting. However, the 
break momentum p^ T 2 is required to be about 45 GeV, 
which is significantly smaller than that of the nuclei. The 
difference of the spectral indices below and above p^ r 2 
is about 0.3. As a comparison, such a value is measured 
to be ~ 0.2 for protons and ~ 0.3 for Helium [lj|. Note, 
for the fit of proton spectrum in a wider energy range, 
as shown in Fig. [I] the spectral difference is only about 
0.1. 

It is also possible that the spectral hardening is not 



a break but a smooth hardening instead, as shown in 
many models [2l|, [25( . We may use a log-parabolic func- 
tion to approximate the smooth hardening of the electron 
spectrum. Fig. 2] presents the results of the fit with log- 
parabolic shape of the primary electron spectrum. We 
find the fit is also improved, with the minimum x 2 value 
slightly larger than that with p^ r 2 free. The fitting pa- 
rameters are given in Table [I] 

From above we see that including a spectral hardening 
of the primary electron spectrum, both the PAMELA, 
AMS-02 and Fermi data can be well fitted with a single 
component of the extra sources. It is a natural expec- 
tation that there is a hardening in the primary electron 
spectrum, given the observed hardening of the CR nu- 
clei. However, the position of the break might be dif- 
ferent from that of nuclei. It is a problem needs to be 
further understood theoretically. If the hardening of the 
primary electron spectrum can be confirmed, it would 
be important to understand the origin and acceleration 
of the Galactic CRs. Since AMS-02 could measure the 
pure electron spectrum with much higher precision that 
PAMELA, we give the expected pure electron spectra in 
Fig. [5] for the above three cases of the hardening. The 
future AMS-02 data may test the existence and detailed 
shape of the primary electron spectrum. 

Finally we discuss the DM annihilation as the sources 
of the e . The annihilation final states are assumed to 
be /i + /i~ and t + t~ . The primary electron spectrum is 
parameterized with Eq. (1), and p^ r 2 is allowed to be free 
in the fit. The fitting results are shown in Figs. [5] and 
for /i + /i~ and t + t~ final states respectively. The fitting 
parameters are compiled in Table [TTJ It is shown that 
the DM models can give comparable fittings to the data 
compared with pulsars. The break momentum of the 
primary electrons, p^ 2 is also similar with that derived 
in the pulsar scenario, and is smaller than p^ r 2 « 230 
GeV. 

The lcr and 2a favored regions on the m x — {uv) pa- 
rameter plane are given in Fig. [51 For /i + /i~ channel DM 
with mass 1 — 2 TeV is favored, while for t + t~ channel 
the mass is obtained to be 2 — 4 TeV. The boost factor of 
the annihilation cross section compared with the natural 
value to give the proper relic density is about hundred to 
thousand. Such results do not differ much from the ones 
obtained through fitting the PAMELA positron fraction 
and the Fermi/HESS e ± data [H. 

The exclusion limits on the DM annihilation into /i + /i~ 
and t + t~~ pairs by 7-rays from the Galactic center (thin 
lines, [50]) and the dwarf galaxies (thick lines, [5l| ) are 
also plotted in Fig. [5) The results show that for the 
t + t~ channel the Fermi 7-rays always give very strong 
constraints on the annihilation cross section. The con- 
straints for the yU + /i _ channel is weaker. The Galactic 
center 7-rays tend to exclude the parameter space to ex- 
plain the e ± excesses. However it may suffer from the 
uncertainties of the density profile of DM in the halo 
center. The more robust limits from the dwarf galaxies 
can not exclude the favored parameter region. 
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FIG. 2: The positron fraction (left) and electron spectra (right) for the background together with a pulsar-like component of 
the exotic e . The high energy hardening of the primary electron spectrum is approximated with a broken power-law and the 
break momentum p£r 2 is fixed to be ~ 230 GeV. References of the data: positron fraction — AMS01 fil. H EAT94+95 J47J ]. 
HEATOO 0, PAMELA 0, AMS-02 Q]; electron — PAMELA [jj, ATIC g<J, HESS [2|I|, Fermi-LAT [n}. 
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FIG. 3: Same as Fig. [2] but the high energy break of the electrons is relaxed in the fit. 
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FIG. 4: Same as Fig. [3] but for a log-parabolic approximation of the spectral hardening of the primary electron spectrum. 



IV. CONCLUSION AND DISCUSSION 

The study of the highly precise data of positron frac- 
tion in CRs reported by AMS-02, as well as the pure 



electron spectrum measured by PAMELA and the total 
e spectra measured by Fermi and HESS, shows that it 
is difficult to use a single component of the extra sources 
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FIG. 5: Expected total fluxes of the pure electrons for the 
three fits corresponding to Figs. [2] - [4j 



to explain the e^ excesses (Paper I, |14j). In this work 
we show that an additional break of the primary elec- 
tron spectrum can improve the fit significantly. The best 
fitting break momentum is about 40 — 60 GeV and the 
spectral difference 72—71 is ~ 0.3 — 0.4. As a comparison, 
the break momentum of protons is about 230 GeV, and 
the spectral difference is ~ 0.1. The hardening behavior 
of the electron spectrum is different from that of nuclei, 
which makes the understanding of the fine structures of 
the CR spectra more difficult. 

In the presence of a hardening of the primary electron 
spectrum, both the pulsar and DM scenarios can give 
comparable fit to the data. However, the DM scenario 
are strongly constrained by the 7-rays, especially for the 
tauon final state. We would like to point out that it 



will be equivalent to take the harder part of the electron 
spectrum and to drop the constraints from the PAMELA 
electron data. In such ways both the AMS-02 positron 
fraction and Fermi total e spectrum can be fitted simul- 
taneously. 

The AMS-02 will measure the electron spectrum with 
high precision in the near future. Whether there is a 
hardening in the electron spectrum or a lower e total 
spectrum than Fermi will soon be answered by AMS-02. 

Note — During the final stage of this work, similar 
idea to improve the fit has been raised in [16| . 



Appendix A: Results of the background positrons 
and electrons 



For the convenience of use we tabulate the fluxes of the 
background positrons and electrons calculated with the 
best fitting parameters of the pulsar models in Table Mil 
For DM models the results have little difference. Note for 
background positrons an additional factor c e + as given in 
Table U needs to be multiplied. Here the local interstellar 
fluxes are given. If one wants to better reproduce the 
low energy electron spectrum, the solar modulation with 
modulation potential given in Table U is necessary. 
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FIG. 6: Same as Fig. [3] but for the DM annihilation into a pair of muons as the extra source of the positrons and electron. 
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FIG. 7: Same as Fig. [3] but for the DM annihilation into a pair of tauons as the extra source of the positrons and electron. 



TABLE I: 


Fitting results of pulsar-like model for different primary electron spectra 




Pbr,1 


= 230 GeV 


Pbr,2 free 


log-parabolic 


Xmin/dof 


236.6/150 


159.0/149 


177.1/150 


parameters 


best 


mean 


best mean 


best mean 


log(A e a ) 


-8.973 


-8.974 ± 0.005 


-8.972 -8.978 ± 0.005 


-8.982 -8.980 ± 0.004 


7o 


1.501 


1.515 ±0.014 


1.526 1.551 ±0.035 


1.503 1.535 ±0.027 


log( P g r , 1 /MeV) 


3.580 


3.596 ± 0.021 


3.589 3.608 ±0.021 


3.620 3.638 ± 0.020 


7i 


2.678 


2.682 ±0.011 


2.773 2.780 ±0.015 


3.178 3.184 ±0.042 


log(pL, 2 /MeV) 


— 


— 


4.776 4.730 ± 0.066 


— — 


72 


2.359 


2.371 ± 0.033 


2.503 2.512 ±0.027 


0.091 0.091 ±0.007 


log^psr") 


-24.913 


-25.024 ± 0.142 


-25.807 -25.753 ±0.118 


-25.802 -25.730 ±0.062 


a 


1.904 


1.884 ± 0.028 


1.742 1.754 ± 0.023 


1.737 1.752 ±0.014 


log(p c /MeV) 


6.243 


6.306 ± 0.066 


6.280 6.301 ± 0.069 


6.130 6.201 ±0.084 


c e + 


1.400 


1.442 ± 0.074 


1.894 1.866 ± 0.070 


1.730 1.753 ±0.051 


<P/MV 


555 


550 ± 28 


650 645 ± 31 


595 596 ± 24 



"Normalization at 25 GeV in unit of cm 



^MeV" 



^Normalization at 1 MeV in unit of cm 3 s 1 MeV 
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TABLE II: Fitting 


results of DM annihilation model 






n + n~ 




r+r- 


Xmin/dof 


188.1/150 


161.3/150 


parameters 


best 


mean 


best 


mean 


log(A e °) 


-8.971 


-8.969 ± 0.006 


-8.981 


-8.979 ± 0.008 


7o 


1.745 


1.767 ±0.037 


1.725 


1.697 ± 0.049 


log(p br ,i/MeV) 


3.675 
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FIG. 8: la and 2<r parameter regions on the m x — (av) plane 
for the DM annihilation scenario. The lines show the 95% 
upper limit of Fermi 7-ray observations of the Galactic center 
(thin lines, with different normalization of the local density 
corrected, [5Cj]) and dwarf galaxies (thick lines, |5l[) for fi + fi~ 
(black solid) and t + t~ (blue dashed-dotted) channels respec- 
tively. 
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TABLE III: 


Local interstellar fluxes of the background positrons and electrons 


5(GeV) 




5 2 5(GcV m" 


2 s- 1 sr- 1 ) 






e+ 


e - (Pb r ,2 = 230GeV; 


e"(Pbr,2 free) 


3 (log-parabolic) 


1.00005 - 01 


1.00855 + 01 


6.37285 + 01 


8.55445 + 01 


7.00795 + 01 


1.26495 - 01 


1.33435 + 01 


8.17495 + 01 


1.09475 + 02 


8.99765 + 01 


1.59995 - 01 


1.6996B + 01 


1.02635 + 02 


1.36785 + 02 


1.12975 + 02 


2.02365 - 01 


2.09395 + 01 


1.26495 + 02 


1.68705 + 02 


1.39415 + 02 


2.55955 - 01 


2.48805 + 01 


1.52875 + 02 


2.03125 + 02 


1.68095 + 02 


3.23755 - 01 


2.86815 + 01 


1.80265 + 02 


2.38885 + 02 


1.98505 + 02 


4.09495 - 01 


3.20945 + 01 


2.07765 + 02 


2.73345 + 02 


2.27535 + 02 


5.17955 - 01 


3.46045 + 01 


2.30515 + 02 


3.02465 + 02 


2.52845 + 02 


6.55135 - 01 


3.57815 + 01 


2.47435 + 02 


3.22885 + 02 


2.71245 + 02 


8.28645 - 01 


3.51825 + 01 


2.54255 + 02 


3.30515 + 02 


2.78865 + 02 


1M81E + 00 


3.21145 + 01 


2.48215 + 02 


3.19665 + 02 


2.71465 + 02 


1.32575 + 00 


2.72855 + 01 


2.31215 + 02 


2.94505 + 02 


2.51575 + 02 


1.67685 + 00 


2.12385 + 01 


2.04835 + 02 


2.57755 + 02 


2.22845 + 02 


2.12105 + 00 


1.49435 + 01 


1.73915 + 02 


2.15935 + 02 


1.88625 + 02 


2.68275 + 00 


9.62365 + 00 


1.43205 + 02 


1.74835 + 02 


1.55405 + 02 


3.39325 + 00 


5.77135 + 00 


1.14735 + 02 


1.37575 + 02 


1.25655 + 02 


4.29195 + 00 


3.31595 + 00 


8.56465 + 01 


1.02165 + 02 


9.61945 + 01 


5.42875 + 00 


1.93675 + 00 


6.03265 + 01 


7.05635 + 01 


6.69455 + 01 


6.86655 + 00 


1.18055 + 00 


4.20105 + 01 


4.80535 + 01 


4.54965 + 01 


8.68515 + 00 


7.59885 - 01 


2.95175 + 01 


3.29755 + 01 


3.11255 + 01 


1.09855 + 01 


5.12305-01 


2.09825 + 01 


2.28365 + 01 


2.15645 + 01 


1.38955 + 01 


3.58595 - 01 


1.51575 + 01 


1.61485 + 01 


1.53595 + 01 


1.75755 + 01 


2.54705 - 01 


1.09645 + 01 


1.15225 + 01 


1.09635 + 01 


2.22305 + 01 


1.83155-01 


7.98315 + 00 


8.21815 + 00 


7.92445 + 00 


2.81185 + 01 


1.33795 - 01 


5.84915 + 00 


5.92875 + 00 


5.79155 + 00 


3.55655 + 01 


9.75975 - 02 


4.27265 + 00 


4.28605 + 00 


4.22885 + 00 


4.49845 + 01 


7.33975 - 02 


3.19855 + 00 


3.18495 + 00 


3.18175 + 00 


5.68995 + 01 


5.35495 - 02 


2.30585 + 00 


2.32425 + 00 


2.31635 + 00 


7.19695 + 01 


3.98775 - 02 


1.69815 + 00 


1.76195 + 00 


1.72335 + 00 


9.10305 + 01 


2.97335 - 02 


1.24955 + 00 


1.34165 + 00 


1.27945 + 00 


1.15145 + 02 


2.22655 - 02 


9.23865 - 01 


1.02015 + 00 


9.60545 - 01 


1.45635 + 02 


1.65585 - 02 


6.81495 - 01 


7.71405 - 01 


7.16735-01 


1.84215 + 02 


1.23495 - 02 


5.10925-01 


5.85515-01 


5.39435 - 01 


2.33005 + 02 


9.22905 - 03 


3.87795 - 01 


4.43645 - 01 


4.06985 - 01 


2.94715 + 02 


6.84055 - 03 


3.01615-01 


3.33965 - 01 


3.06645 - 01 


3.72765 + 02 


5.08845 - 03 


2.35015 - 01 


2.52155-01 


2.32255 - 01 


4.71495 + 02 


3.76415 - 03 


1.82535 - 01 


1.89315-01 


1.76015 - 01 


5.96365 + 02 


2.76415 - 03 


1.41015-01 


1.41875 - 01 


1.33905 - 01 


7.54315 + 02 


2.03085 - 03 


1.09545 - 01 


1.06575 - 01 
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9.54105 + 02 


1.47385 - 03 
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7.94855 - 02 
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3.03915 - 02 
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1.35415 - 02 
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8.00385 - 05 
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1.06895 - 02 


7.90605 + 03 


5.04135 - 05 


8.37325 - 03 


5.83845 - 03 
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